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Chern-Simons gauge theory, since its inception as a topological quantum field theory, has proved 
to be a rich source of understanding for knot invariants. In this work the theory is used to explore 
the definition of the expectation value of a network of Wilson lines - an embedded graph invariant. 
Using a generalization of the variational method, lowest-order results for invariants for graphs of 
arbitrary valence and general vertex tangent space structure are derived. Gauge invariant operators 
are introduced. Higher order results are found. The method used here provides a Vassiliev-type 
definition of graph invariants which depend on both the embedding of the graph and the group 
structure of the gauge theory. It is found that one need not frame individual vertices. However 
without a global projection of the graph there is an ambiguity in the relation of the decomposition 
of distinct vertices. It is suggested that framing may be seen as arising from this ambiguity - as a 
way of relating frames at distinct vertices. 
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I. INTRODUCTION 



o 
o 

Topological quantum field theory, developed by Atiyah |jj and Witten ||, is rooted in the desire to construct 
ON i a framework independent of background structure. Chern-Simons gauge theory, being diffeomorphism invariant, 
provides an ideal test case. Witten |3| made the remarkable step of relating the vacuum expectation values of Wilson 
loops in Chern-Simons theory to the Jones polynomial [Q. This result was expanded to encompass more general 
O 1 groups and representations and observables based on projected graphs by Witten || and Martin ||. Parallel to this 
work, techniques of non-perturbative, background independent quantization were developed |Q - ]l5| |. Applied to 
gravity these techniques have led to a thorough understanding of quantum geometry Jl6| (Tt]] . It was discovered that 
in order to describe quantum geometry it is necessary to consider not only states based on projected Wilson lines with 
• ■ *h , intersections but also to consider the full three dimensional spatial structure of vertices; the tangent space structure 
/\ ' of embedded graphs is required (l6) , JlTj . 

This paper, using the variational technique introduced in (Ts| , generalizes the early results of Witten |j| and Martin 
H to embedded graphs with vertices of arbitrary valence and general tangent space structure. It is seen that the new 
invariants introduced here contain a dependence on the relative orientations of edge tangents. Further, it is seen in 
detail how the first order variation exponentiates to the full, non-perturbative results expressed in terms of Temperley- 
Lieb recoupling theory jUl (in certain projections). With these techniques it is easily seen that the balanced networks 
of Barrett and Crane |20|| are easily seen to arise from the first order formalism. Finally, the variation calculations, 
accounting for the tangent space structure at vertices, strongly suggests that framing, originally introduced as as way 
to partially restore diffeomorphism invariance, is gauge. 

This study comes from the confluence of the two approaches: the calculation of vacuum expectation values in Chern- 
Simons theory such is in Refs. ||, [||, |2l), 22 1, and |2^] and the development of background independent quantization 



techniques - jl5). The relation between this fields is actually tighter. More than a decade ago Kodama noticed 
that the Chern-Simons form is a formal solution to all the constraints of the vacuum theory with a cosmological 
constant [^ij . As it is also natural to consider loops in this context as well, expectation values of Wilson loops in 
Chern-Simons theory have a dual role as "loop" transforms of this Kodama state in the connection representation of 
quantum gravity. 
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As we know that expectation values of loops in Chern-Simons theory require framing || the question arises as to 
what is the freedom associated with a framed vertex. The variational method provides a key to explore the definition 
of such singular graphs pj| . As shown in this paper the character of the invariant changes since an interplay of group 
and manifold structures determine its value. In fact, one recovers the Temperley-Lieb recoupling theory of Kauffman 
and Lins ]l9| tempered by a dependence on the embedding of the graph. 

Motivation for this work also comes from a desire to have a more complete understanding of the loop representation 
for the Kodama state. Realizing that the loops had to be framed, a new representation of loop algebra of quantum 
gravity was introduced in Ref. |2q| . In the framed loop representation, products of operators must also be defined. 
Two questions arise immediately: Is the product uniquely defined? If the product of operators is not unique, what is 
the freedom due to framing at a vertex? Thus it is natural to investigate the role of framing play at a vertex. From the 
perspective of the framed loop representation or "^-quantum gravity," it is interesting to find the vacuum expectation 
value of two loops (or generically graphs) which intersect. One of the goals of this work is to determine whether 
additional structure is required to specify this expectation value. This in turn should be reflected in the basic algebra 
of the operators in g-quantum gravity. 

The main tools used in this study are based on spin networks. Originally introduced by Penrose as a method to 
solve the Four Color Theorem, and then as a combinatorial basis for spacetimc p7[ , spin networks have found a 
new life as a basis for 3-geometry states.^] In this new role spin networks are labeled graphs embedded in the three 
dimensional spatial slice. The spin networks (or "spin nets") resolve a long standing problem of the over completeness 
of the loop representation ||. (A state space built simply from loops is subject to a number of identities called the 
Mandelstam identities.) This new basis solves these identities and comes with a bonus. The spin net basis is the 
eigenspace for geometric operators such as area 16 and volume ]T?j] . The study of the kinematic states of quantum 



gravity using this basis I call, for the purpose of this paper, spin network geometry. 

As Chern-Simons theory is a diffeomorphism invariant theory, it is not surprising that the vacuum expectation 
values of loops are functions of diffeomorphism invariant classes of knots and links (up to the subtleties of framing) . 
However, the presence of singular knots or, in general, vertices significantly complicates the picture. One pleasing 
possibility is to define such graph invariants in terms of non-intersecting links. For instance, it may make sense to 
view intersecting knots as the limit points of non-intersecting knots; the freedom in taking the limit would encode 
precisely the freedom in the constructing intersections. One way to explore these ideas is through the study of 
Vassiliev invariants, in which one defines singular knots by associating to them an invariant defined by a difference 
of the possible limits of non-intersecting knots. For instance, the Vassiliev invariant associated a simple intersection 
is simply the difference between the over- and under-crossing decompositions, i.e. ( X ) = ( X ) — ( X )■ For more 
than one vertex, this may be continued iteratively [^l] . More generally, Vassiliev invariants allows one to associate to 
every knot an infinite sequence of rational numbers. This sequence divides up into finite sub-sequences which form 
vector spaces and, when indexed by an integer n, are Vassiliev invariants of order n. These invariants have a number 
of nice properties. They are related to knot polynomials: By replacing the variable of a polynomial invariant by e x , 
the coefficient of order n is a Vassiliev invariant of order n. Vassiliev invariants of finite type vanish at order i for 
knots with more singular points, Vi(K J ) = if j > i where j is the number of singular points. These form an algebra, 
so that two invariants of finite type, say i and j, yield an invariant of the product ij. Finally and spectacularly, it is 
conjectured that a complete set of knot invariants may be built from Vassiliev invariants @). The variational method 
used here is naturally associated to differences in non-intersecting knot invariants. It turns out that the lowest order 
results are Vassiliev invariants of finite type. 

Recently Gambini, Griego, and Pullin, building on earlier work by Alvarez and Labastida [p2|], have proposed that 
Vassiliev invariants, generalized to include spin net states, are solutions to the Hamiltonian constraint of quantum 
gravity pjJ]. This construction is based on the idea that the framing dependence may be collected into an overall 
factor in the expectation value of Wilson loops. They find that the loop derivative operator is well defined when 
acting on Vassiliev invariants, suggesting that it may be possible perform dynamical calculations entirely in the spin 
net representation. 

In Section || I give a review of the variational technique. This, as in the work of Labastida and Perez M] , associates 
to every vertex a gauge invariant operator by analyzing the relations between families of (non-diffcomorphically related) 
graphs. Such relations are generated by taking the variation of the vacuum expectation value of loops. These loops 
depend on a parameter which interpolates between intersecting and non-intersecting loops. The variation allows one 



to construct an operator for vertices. Section III is devoted to applying this technique to graph invariants. As these 



operators a quite similar to the operators of spin net geometry, the same recoupling techniques may be used to find 



1 The structure for three dimenionsal SU(N) Chern- Simons theories was also outlined in Refs. H and 
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the expectation values. In Section [V I present the results of these calculations. In some cases spin networks prove to 
be eigenspaces of the the operators. In the final section I offer some concluding remarks. 



II. TECHNIQUES OF THE VARIATIONAL CALCULATION 



In this section I give a self-contained review of the variational technique. Though the technique is widely known, 
the presentation serves to fix notation and to emphasize the embedding dependent elements of the calculation. In 
this way it becomes clear how one may generalize the technique to a projection-independent one. The method is an 
extension of those in Refs. f25|| , pjj , and [|T). It relies on several key properties of Chern-Simons gauge theory and 
the definition of spin nets which I mention before performing the variation. In the next section I collect the results 
and define gauge invariant operators for arbitrary graphs. 

Let us consider Chern-Simons gauge theory on a smooth three manifold £ without boundary 

S[A] = f Tr [A A dA + \A A A A A] (1) 

in which A a = A % a T l is the Lie algebra- valued connection (a, 6, c, ... are abstract spatial indices). The gauge group is, 
for this and the next section, taken to be a compact, semi-simple group G. The trace is taken in the fundamental 
representation. The generators T l , i = 1,2,... dim G, are normalized so that Tr [T l T J ] = —^S^. From the perspective 
of canonical quantum gravity, this action has another role as a state in the connection representation. 

Kodama noted that there exists a state which formally satisfies all the constraints of vacuum canonical quantum 
gravity with cosmological constant |24| 

$[A] = exp{iS[A}) 

where S[A] is the Chern-Simons action of Eq. (|l|) with k = 247r/A (A being the cosmological constant). In this 
perspective the manifold £ is the spatial slice in the (3 + l)-decomposition. 

Given the gauge and diffeomorphism invariances of the theory, Wilson loops are natural observables. The vacuum 
expectation value of a (knotted) loop K is related to the Jones polynomial via (3| 

(W K {A}) = q-i w «J K (q). 

The Jones polynomial Jk(q) is the ambient isotopy polynomial invariant of q — exp(7rz/£;). (Omitting, for the present, 
the non-perturbative shift in k. For more details see the appendix.) The writhe, u>k , is given in terms of the sum 
of crossings in a knot diagram of K. In fact, the right hand side is defined only for a projected knot in blackboard 
framing]^] The goal here is to generalize the observable in two ways. First, the variation technique reveals that it is 
possible to remove the projection dependence. Second, as is clear from both gauge theory and spin net geometry, it 
is wise to include observables based on graphs or spin nets. 

Spin nets are embedded graphs labeled (or colored) with integers which represent the number of lines running along 
an edge and, equivalently, identify the irreducible representations carried by the holonomy. Every vertex contains 
a combination of Clebsch-Gordon coefficients which is called an intertwiner. For the purposes of this paper, a spin 
network TV consists of the triple (G; i, n) of an oriented graph, intertwiners, and edge labels. The corresponding spin 
net state S is defined as 

(A|G;i,n)= H i v o (g) U [ne) [A}. 

tiGv(G) eSe(G) 



2 Framing conventions are easy to understand with White's theorem. This states that this self-linking number is the sum of 
the writhe and the twist fj(i) |. Pictorially, self-linking is the winding number of the frame around the base loop; twist records 
the number of sides of the ribbon one sees in a projection (Mobius bands are ruled out); and writhe is given by the number of 
curls in a line. There are a number of framing conventions which fix writhe and/or twist including blackboard and standard 
framing. Blackboard framing, in which the frame is always normal to the knot in the plane of the blackboard, sets the twist 
to zero. By White's theorem, the contribution to the self-linking number of a twist may be expressed as an equal number 
of curls. Standard framing requires the self-linking to vanish in any projection. It is naturally selected since it removes the 
explicit projection dependence. However, this choice only exists for certain manifolds including S 3 H. 
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These states are gauge invariant as the intertwiners are invariant tensors on the group. In more picturesque language, 
the state is gauge invariant because the intertwiners connect all the lines at each vertex. 

To investigate the definition of the vertices, it will be convenient to analyze sub-spin nets. One may view the 
following analysis as cutting out a ball around a vertex, operating on it and then reinserting the result back into 
the graph. The operations are local so it is convenient to keep the intersections of the spin net and the vertex with 
the 2-sphere "fixed." Interestingly, the variational method does not allow one to view these manipulations as solely 
operations on an abstract graph; one must also include tangent space information at the vertex. In this sense the 
calculation is tied to the manifold structure. The sub-spin nets I consider here consist of a vertex, incident edges, an 
intertwiner and the edge labels. This will be denoted as v = (v, e v ;i v , n) G (G; i, n). 

To analyze subgraphs will be necessary to work with Wilson lines or holonomies. I will take a path a to be an 
oriented, piecewise smooth map from the interval I — [0, 1] into S. The composition of two paths will be denoted 
with o as in ei o e2- I will take paths to be non self-intersecting.R Associated to a path one has an holonomy 



U a [A] = V exp f dta a A a 
Jo 



(«(*))• 



Paths with boundary points identified are called loops. Wilson loops are simply the trace of holonomies based on 
loops. For a spin network A/", the vacuum expectation value is defined by the functional integral 

(W M [A]) = [ [DA]e iS ^WMA}. 



Two identities work hand in hand to make the variational calculation possible. The first shows that it costs curvature 
to differentiate the action with respect to the connection 

— ^-cS[A] = ^e abc Fl c {x) (2) 
oA l a (x) STT 

where F a b is the curvature of the connection A a . The second identity concerns the variation of paths. I take families 
of paths (frequently pairs of edges) of a spin network to be parameterized by a continuous parameter usually denoted 
u. For some value of this parameter it, the path intersects other edges. 

Under the variation of a holonomy with respect to this parameter one discovers a magical property. The variation 
costs curvature. The change of an edge, e u labeled by n, of spin network state S is given by 



d 

du 



U eu = dte a (t)e b (t)F l ab (e(t)) U e {Q,t)T( n) U e {tA) (3) 



where e a denotes the derivative of the edge with respect to the parameter u. The original spin net is recovered for 
one value of the deformation parameter uq when e u \ u=Uo . (I will not explicitly show the the dependence of the edge 
e on u when it is clear from context.) 

There are two distinct forms of variations of Eq. (||) |n| . In the first form, used for the decomposition of vertices, the 
paths depend on the parameter u which determines the "location" of the path relative to another, e.g. as in an "over-" 
and "under-" crossing. For this reason, I call u the "decomposition parameter." Though the map e„ is continuous 
in the manifold (like a homotopy transformation) , given the topological nature of the theory, one expects that this 
variation of (Wj\f) to be discontinuous. In fact one finds that the variation of the expectation value is distributional. In 
the second form, in which the loops are parameterized such that the affect is that of a diffeomorphism, the expec tation 



value is naively expected to vanish. As the theory requires framing, this is not always the case (as in Sec. IV C ). The 
difference of these two variations may be simply expressed asu = u(t) for the first and u = u(x) for the second. 
The final piece of the variational calculation is the property of the holonomy 



5 f 

— — U a [A} = / dta a (t)S 3 (x,a(t)) U a {0 t t)Tf n) U a {t,l) 

0A a\ x ) JO L 



(4) 



in which the path a carries the label n. The variation inserts a group generator in the n/2 representation into the 
holonomy at the point x. 



3 There is no loss of generality; if, for instance, a loop had a single self-intersection then, expressed as a graph, the loop would 
have two edges. 
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As a first example of the variational calculation, consider the four valent vertex, or double point, shown in Fig. (Q). 
This is part of a larger spin net Af so that the vertex v = (v, e; i v , n) 6 N. Only the edges e\ and are parameterized 
by u. Further, suppose that this dependence is coordinated so that when u > (u < 0) the intersection is resolved 
into an over (under) crossing as seen when the vertex is projected along the e a direction. When u = the edges 
intersect, forming the double point. 



e-, 

(a.) 




e 3 
(c.) 



FIG. 1. The four valent vertex considered here has four edges ei,e2,es, and a labeled by integers a and b. It is resolved 
into two crossing diagrams as shown in (b.) for u > 0, (Ue 1 oe 3 Ue 2 oe i ) + , and (c.) for u < 0, {U ei oe 3 Ue 2 oe i ) _■ The intersection 
(a.) occurs when u = 0. The u-derivative of the edges, which determines the meaning of the over and under crossings is shown 
in (a.). 

Using the identities of Eqs. and (||) and integrating by parts, the variation of the expectation value ([/vL4]) 
may be expressed as 



du 



(U V [A}) 



[DA]e 



iS[A] 



I dt6t(t)e\{t)e abc J —^ (u ei (0,t)T{ a) U ei (t,l) U e2 U e3 U ei ) 



+ / dte%(t)e b 3 (t)e abc 



U ei U e2 U e3 (0,t)T(a)Ue 3 (tA)Ue 



(5) 



in which U e represents the whole holonomy along the edge U e (0, 1). The variation produces a first order (in 1/fc) 
result with a generator inserted into the network and differentiation with respect to the connection along the edges 
ei and e^. In carrying out this variation, one clearly must make the critical assumption that the derivative commutes 
with the integration over the connection. It is also well to note that the integration over the loop parameter t is only 
over the loop space where e u (t) is non- vanishing. 

The differentiation in Eq. (g) acts on all the holonomies of the incident edges. Thus, with the identity (^) , Eq. (|5|) 
becomes (rti = 713 = a and 712 = 714 = b) 



— (in - — 

du { v) k 



DA]e iS\A] J- / dt^(t)e$(t) (0,^)^(4,1) 



4 1 4 

[ dseJ( S M 3 ( ei (t),e fc ( S ))C/ efe (0, S )T/ nji) ?7 efe ( S! l) f[ U, 
fc=i,fc#i Jo i^j,k 



(6) 



As complicated as this might appear, the structure is rather simple. One gauge generator is inserted on the edges 
which depend on the parameter u. The other generator is inserted in all other edge pairs, generating the second sum. 
Further, the delta-functions deflate these terms to terms specified by the condition ej(t) = efc(s). For example, the 
terms 



f dtel{t)e\{t)e abc [ ds e c 3 (s)S 3 ( ei (t),e 3 (s)) ( U ei (0, t)T{ ni) U ei (t, 1) U e3 (0, s)Tt n3) U e3 (s, 1) J] U ei 

Jo Jo \ 1^1,3 

abc [ d S ef( S )<5 3 (e 3 (4), ei ( S )) (^(0,5)^^(5,1)^3(0,^)^(4,1) [] U e . 



(7) 



+ / dtei(t)e b 3 (t)e a 



only differ by a sign when the condition is satisfied and thus cancel. 

There are two classes of solutions to the above condition. For a single edge, j = k, there is a one-dimensional 
solution s = t. These singular terms require special care. Us ually the line is split in two to give a loop and its frame. 
I will postpone discussion of this type of term until Section IV C. The second class of solution, for j ^ k, lives only 
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at a vertex. For instance, the first term of Eq. (o) may be expressed as 



47TI 
~k 



dte\{t)e\{t)t ahc / dse^S* ( ei (t),e 2 {s)) (U ei {0,t)T{ a) U ei (t, l)U e2 (0, s)T{ b) U £2 {s,l)U e3 U e4 



(8) 



The 5-function reminds us that the variation changes diffeomorphism equivalence class. The singular nature of this 
term is contained in the "volume" term 

dt J dse abc e a 1 (t)e b 1 (t)e c 2 ( S )6 3 (e 1 (t),e 2 ( S )). (9) 

Nevertheless this result is useful. One may rewrite the delta-function as 

1 



(ei(t),e 2 ( a )) 



|A(i, 2) (0, 0,1)| 



5(u)5(t)5{s - 1) 



in which 



A ( j,fc)(u,f, s) = e afcc e"(w)e5(t)e£(s). 



for edges ej and . Upon integrating over an interval around u = the term of Eq. rtq) may be simply expressed as 



+c 



T{a)U ei 



in which 



k(j, ft) 



A(j, fc )(uo,0,l) 
|A( jjfe )(u o ,0,l) 



(10) 



(11) 



The sign factor k) is only defined at the intersection when u — uq. I also use the convention that dx 5(x) = ^. 
It is clear from the structure of the operator of Eq. jl0| ) that the affect of the variation is simply to act with a 
left or right invariant vector field on the incident edges (indicated with square brackets above). The handedness is 
determined by the orientation of the edge. It is clear that both group structure and tangent space structure determine 
the variation of the invariant. Thus, for this type of solution, the result may be expressed as insertion of generators 
at the vertex times a sign. In the next section I present the form of the variation for this vertex and more general 
vertices. In each case they correspond to lowest order Vassiliev invariants. 



III. GAUGE INVARIANT OPERATORS FOR ARBITRARY VERTICES 

In Ref. gauge invariant operators were constructed for knots with planar double points. From the perspective 
of spin net geometry, in which states of the theory based on planar vertices have vanishing volume expectation values 
p7[ , it is clear that one would like to generalize the construction to non-planar higher valence vertices. This section 
contains an analysis for these cases. The next section contains an evaluation of some of these operators on spin 
network states using SU(2) group structure. 

As an example of the gauge invariant operators for vertices, consider the four vertex calculation presented in the 
last section. The terms on the the right hand side of Eq. (^) , depending on the nature of the volume term, fall into 
two classes. Here I will only study the terms between distinct edges, postponing the "self-interaction terms" until 
Section IV. C. As noted in the last section the two terms with the edge pair e\ and e% cancel, due to the opposite 
signs. Thus, Eq. ffl) reduces to 



du (U v ) — (U ei oe3 Ue 2 oe4 )_|_ \U ei oe3 ^e 2 oe 4 }_ 

= — 2) + k(1, 4) + k(3, 2) + K (3, 4)) 



(12) 



It is understood that the holonomies are over complete edges and the brackets indicate the composition of the 
holonomics. The sign is determined by the decomposition parameter and the tangents of the incident edges. For such 



G 



a simple vertex such as this 4- valent one, all the terms of Eq. (g) combine into one. For higher order vertices this is 
not the case. Clearly the tangent space structure and the choice of the vector field e u (t) determine the decomposition; 
the sign terms collect as overall factors on each term. In this 4-valent vertex, the four terms may take nine possible 
values which are realized by different configurations of the incident edges. For instance, the vertex shown in Fig. (Q) 
has an overall factor of 4. Similarly in the planar case where the edge pairs have the same tangents, all the signs are 
identical and this simply reduces to one term |2l|j 



(±1) 



47TI 



Ue 3 T{ a )U ei 



(13) 



The overall sign is determined by the direction of the vector field associated with the parameter u. Since the sign also 
changes on the left hand side of Eq. ( |l2| ) , there is no additional freedom in determining the planar 4- vertex. 

As this last point will become more important later, it is worth expanding on. The vector field e u determines the 
deformation of the vertex. One may lift one edge pair using such a vector field as is done here. One may also "slide" 
a vertex along an edge by taking a vector at the vertex e„(0) which is parallel with one of the tangents at the vertex. 
In the case of a simple 4- vertex such a slide is redundant since the lines are reparameterization invariant; the spin net 
state does not change.^] However, it illustrates the edge of the region of equivalent decomposition vectors. For the 
4- vertex of Fig. (Q), all vectors e a | M= o in the tangent space volume designated by the triple e\, 64, e^ 1 yield equivalent 
decompositions in that they are diffeomorphic |32j] . The compliment also yields the same decomposition. 

For graphs with more than one vertex, this procedure of lifting edges iteratively may be applied to all the vertices 
in the graph f2l|| . The procedure may also be applied to higher valent vertices. In this case, for a 2n- valent vertex 
one uses n — 1 variations of this type. The lowest order in 1/fc of this decomposition is then also of order (n — 1) 
on account of the integration by parts for each variation. It is worth noting that gauge invariant operators of this 
type may only be associated to vertices with even valence. Otherwise one would be left with an open edge and the 
resulting state would no longer be gauge invariant. 






e 5 o e 2 



FIG. 2. The 6 valent vertex v = (v,e;i v ,n) is shown with an intertwiner tree, (b.). The over-crossings in (b.) only indicate 
the nature of the connections in the intertwiner. The dotted line indicates that the diagram inside has no spatial extent. The 
paths ei o and ei o es are parameterized by u and v, respectively, and varied in the calculation. The first term of the variation 
is given in (c). 



Before giving the general case, I will present the derivation for a 6-valent intersection. Consider the vertex shown 
in Fig. (|^) with six edges and an internal tree with labels and 13. This vertex may be decomposed in at least 
three ways. For instance, the vertex may be decomposed first into a four valent vertex created from e2, e^, e 5 , and e§ 
by lifting the edge pair ei o a. The decomposition may be completed by lifting the edge pair es o e2. As with the 
Vassiliev invariants for double points, the intersection may be expressed as a signed sum of over- and under- crossings. 
Performing such a decomposition by first lifting the line e\ o e4 one finds that the eight terms begin with 



-^{U v ) = ^ I </* / < //'„/„.n(/i.-';(/i^(-)'V ! (r J .(/i.^(,i) 



x (i v o U ei (0,t)T{ a) U ei (t, l)U ei U e5 U e2 (0, s)T{ b) U ea (s, l)U es U e A + 



(14) 



in which the holonomies and generators are composed with the intertwiner. To finish the decomposition one can lift 
a second line, here chosen to be e2 o es. This line is parameterized by v and the vector field e v (t) is taken to be in 



4 For higher valence vertices, even planar ones, such a slide changes the valence of the original vertex and creates more vertices. 
In these cases, the variation is non-vanishing. 
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the same direction as the vector field associated with u. After integration in the two parameters u and v the result 
becomes 

XW-/f w-x 



4vrA 2 

TJ 



K 2 ({1, 4}, {2, 5}) (i v o U^T^T^Ue^T^TgUe^Ue, 



(15) 



-«({1,4},{2,5})«({2,5},{3 J 6}) (i v o U ei T{ a) U e4 U^T^U^T^U, 
-«({1,4},{3,6})«({2,5},{3,6}) (i v ° U ei T? a) U e4 U e2 T'U e5 U e3 T^Tp c) U, 



-«({1,4},{3,6})«({2,5},{1,4}) (i v o u^t^u^u^u^u^u,,, 

(a) U <^i U es T {b) T (b) U ^2 Ue 3 Tfa U e6 
(a) Uet U e2 T^ b) U e5 U^T^T^U^ 

in which the signs have been collected with 

j}, {k, I}) = k) + k(i, I) + n(j, k) + I) 

and T^T^> := T l T^ + T^T 1 - symmetrization defined without a numerical constant. 

A few remarks are in order. First, the diagrams on the left hand side are a schematic representation of the 
decomposition. The plane of the projection is given by the vector e u (it) (it is the value of the loop parameter 
at the vertex). Different choices of the vector fields e v (t) and e u (t) yield different decompositions. For a planar 
vertex, there are three independent ways to decompose the vertex corresponding to permutations of the three paths 
from which the vertex is built. The choice of decomposition is made by which edges are lifted and limits of the 
decomposition parameters. (The order of the lifting does not influence the result as the variations commute.) All 
possible decompositions may be so generated. Here, the limits of the v integration are chosen to be less than the 
u integration, i.e. e u > e v . Second, this result is explicitly second order as is expected for a Vassiliev invariant of 
finite type. Third, the tangent space structure and the group structure separate in each term. Fourth, no framing 
was required in the decomposition because the variation did not produce volume terms which had to be regulated. 
Of course, part of this was by fiat since the edge self-linking terms were neglected (these will be discussed in Section 



IV C). Nonetheless, this calculation does indicate that a vertex decomposition does not need such a regularization. 
Fifth, as both the operators and the intertwiners are invariant tensors on the group, the operators are gauge invariant. 

The iterative procedure used on the 6-valent vertex may be carried out on higher valence vertices as well. An even 
valence vertex may be completely decomposed in this manner while odd valence vertices have a minimally trivalent 
decomposition. The general decomposition of a 2n-valence vertex is conveniently expressed if an incident edge i is 
numbered so that its partner (the edge which joins to the first when the vertex is decomposed) is numbered n + i. It 
will also be convenient to label the pairs by the first edge. For instance, the edge &\ has partner e( n+ i) and the pair 
ei o er n +i) is labeled by the index 1. The general form of the operator is clear from the two previous calculations. 
Generators are inserted in all the lifted edge pairs. The other (n — 1) generators may be inserted in any of the 
remaining (n — 1) edge pairs. Since these permutations are all distinct one may index them with one parameter m. 
There are N = (n — 1)™ _1 possibilities. For the 2n-vertex V2„ = (v, e v ; i v , n) one has 

/ .\ (n-l) N n 

D (n - 1) (U V2 J= [ j) E^n^MtfM^+n))- ( 16 ) 

^ ^ m— 1 e— 1 

The map (f>r n \ (m) gives the generator insertions 

P(n e ) V 11 ) ~ 1 (n <1 ) 1 (n <! ) ■ • • 1 (n e ) ' 

Each possible term, indexed by m, contains a symmetrized set of generators. The sign factor is given by the product 

n-i 



K r . 



'[[ K {{I, {I + U)}, {p m {l), (Pm(l) + tl)}) 



in which p m (l) labels the edge pair induced by the permutation m. For example, the last term of Eq. ( |15|) indexed 
by to = 4 would have P4(l) = pi{2) = 3 while all other values of p± vanish. The overall structure of Eq. (|16|) is easy 
to see: All lifted pairs have generators. The sum is over all the possible insertions of the remaining generators while 



8 



the product is over the incident edges. A decomposition of an entire graph would include a second product over all 
vertices. 

The vacuum expectation value of the operators of Eq. ( pLf ) gives an invariant for the singular graph with a 2n- 
valent vertex. This invariant may be expressed as a signed sum of non- intersecting knots or links. Since the result 
is minimally of (n — 1) order, the invariant vanishes for all graphs with an additional crossing (either in the vertex 
under consideration or at another site on the graph). Thus, this provides an instance (in a slightly different setting) 
of the theorem by Birman and Lin which states that the n-th order coefficient of the expansion of a polynomial knot 
invariant is a Vassiliev invariant of order n |28|j . 

For singular knots and links these operators may be given a graphical form. In these diagrams, each component 
of the link is represented by a circle with marked points. These points represent the values of the loop parameter(s) 
where an intersection occurs. The generators are inserted at these points. The contractions of the group index are 
represented as a dashed line. For the knot shown in Fig. (0), one operator, 



Tr U ei T {i T 3) U e2 T k U ei T l U e , T j U e3 T k 



is shown in part (b.). The pair of indices ij at Sj are symmetrized. This is represented in the figure as a pair of 
crossed dashed lines. For multicomponent links the distinct circles will, generically, have dashed lines between them. 
Since each term in the decomposition is a distinct operator, each diagram defines a distinct configuration. Different 
links may be classified according to their diagrams. The diagrams may be seen to label a product of jV„-dimensional 
vector spaces (where N v is as above for a 2n- valence vertex). A link may be identified by its configuration. These 
diagrams are also useful for computing invariants of expectation values of spin net states with a flat connection, such 
as the numerical invariants of Ref . p2fl . 





FIG. 3. The operator associated to the singular knot given in (a.) may be represented with a diagram, e.g. (b.). The 
operator represented is ^Tt - one of four terms. The tangent space dependence has been 



omitted. 



It is intriguing to note that, between vertices of the same graph, there exits no natural relation between the vector 
fields associated to the decomposition parameters at distinct vertices. Of course, in a global projection of the entire 
graph, there is a natural choice of the vectors fields, such out of or into the page. However, in the absence of such 
a projection, the relation must be consistent with an arbitrary projection at each vertex. What one might seek is a 
way to leave the projection choice arbitrary at each vertex and yet still have a consistent system. Such a system is 
reminiscent of a gauge. 



IV. EVALUATING 517(2) VERTEX OPERATORS 



For SU(2), it is relatively simple to express the invariant operators of the last section in terms of spin networks. The 
techniques are similar to those used in the geometric operators of quantum gravity (see jl6| and [0). As in the case of 
these operators it is convenient to use the methods of recoupling theory. In this section I shall present the calculations 
using the diagrammatic methods of Kauffman and Lins Jl9| ). Of course, the recoupling theory is for the group SU(2) 
not for the quantum group SU{2) q . When the spin network basis is an eigenbasis for the gauge invariant operators 
(as they are in some cases), the variation operator may be exponentiated to give the full series. The advantage is 
immediate. The simple relation between Temperley-Lieb recoupling theory and the "classical" or "binor" conventions 
of regular SU(2) recoupling theory, in which ordinary recoupling theory is obtained when A — > —1 |19|, suggests that 
the variation yields the lowest order result of the invariant, "The classical result gives the quantum exponent." 
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To find the action of the SU (2) graph observables it is necessary to introduce a spin net decomposition of the gauge 
generators. This extension of the simple relation 



[T*]JJ IT 



\c 



(17) 



(for SU (2) using binor conventions) for single lines becomes 

a+b 



T 



(a) 



T 



= 7 E (-l) MQ ' b) «c(a,6)- 



c— | a — b\ 



(a, 6, c 




(18) 



in which 



a c (a, b) = - [a(a + 2) + 6(6 + 2) - c(c + 2)] . 



The recoupling quantities 0(a, 6, c) and A„ are given by Eqs. ([34|) and (j3ll) in the Appendix.^] Dotted ovals in these 
diagrams represent the recoupling at the vertex. The identity of Eq. (|18|) may be applied to the gauge invariant 
operators of the last section. I will give results for the 4-valent, 3-valent, and 2-vertices before turning to a brief 
discussion of the general case. 



A. Four valent vertices 



With the above identity in hand, one may compute the action of the operators of the last section on the spin net 
states. I will only evaluate the operator for planar vertices here. The extension of these results to arbitrary four-valent 
vertices is straightforward. 

The variation of the planar four-valent vertex of Eq. (13) with the intertwiner 'X indexed by i is 



= S£(-ir* (o ' 6) acM) A 



9(a, 6, c) 




™ ( _ 1)a ,(a, 6) ^(« + 2) + 6(6 + 2)-»(z + 2) j / 




(19) 



in which the identity (Q) was used in the first line. The recoupling calculation of the first term is simple with the 
use of Eq. ( |35|) while the recoupling for the second terms is more involved, using identities (|36|), (|42]), and (^). It 
is interesting to note that spin nets are an eigenbasis for the variation. This suggest that one ought to exponentiate 
the first order result. The exponentiation is not unique. However, by examining the sum over intertwiners, one can 
match these results with those of Temperley-Lieb recoupling theory. 

When evaluating the operators such as the 4-valent vertex operator of Eq. dl3|), the vertices are labeled by 
intertwiners. If we are to determine the expectation value of the product of spin net states, which have an pointwise 
intersection, the calculation is not defined from the group theory standpoint until some sort of intertwiner is specified. 
When an intertwiner is not specified it is most natural to sum over the possible intertwiners as in Eq. (|45[). The 
variation then yields 



= S £(-l) (a+b -* )/2 Ha + 2) + 6(6 + 2) - i(i + 2)] 



#(ct, 6, i) 



(20) 



(which is identical to the variation without specifying an intertwiner). This agrees with the result from Temperley- 
Lieb recoupling theory in so much as the first order coefficients are equivalent: Using Eq. (fl5h the two crossings may 



5 The norm used in the identity differs from the normalized spin networks of Ref. jL4| by a factor of \J A c 
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be expressed in terms of the intertwiner 



(X)-(X)-E^-^,-)^(g) 

Expanding the A coefficients to first order one finds 

( X ) - ( X ) - (£«,<«, «) ^ ( g) ) + °<i/*"> 

which matches the variational result of Eq. (|2(]). In this comparison one learns that the first order variation captures 
the first order term and the overall sign. Thus, to match the conventions of Temperley-Lieb recoupling theory, one 
must exponentiate only the first order dependence on the labels, not the overall sign. 

When evaluating the vacuum expectation value of a product of operators based on knots which intersect, there is 
not a unique choice of intertwiners for the new vertices. The most natural solution is to sum over possible intertwiners. 
The above calculation shows that a product of operators based on intersecting loops will be a sum of operators. Thus, 
in the perspective of g-quantum gravity, a product of intersecting loop operators will generically give a sum over 
admissible states. It may be best to use a new basis or a new set of operators based on the set of eigenstates of the 
T q operator of Ref. |||. 

In is interesting to note that the variational technique produces only a restricted set of invariants associated with 
signed sums of vertex decompositions. For instance, the relation for single lincs^] 

X =- a *+a-* (X + X ) 

used in Ref. [^6) as the "deformed Mandelstam identity" simply does not appear. This "averaging decomposition" 
clearly is an invariant of a different character than the Vassiliev invariants derived here. However, this calculation 
allows for the possibility that the Vassiliev invariants and the invariants associated to this averaging procedure both 
have an interpretation in terms of the expansion of the polynomial invariant. The expansion of the average invariant 
begins with terms of Oth and 2nd order while the Vassiliev invariant begins at the 1st order. The two sequences could 
live at different orders in the expansion of the polynomial invariant. The average method of vertex decomposition 
could perhaps be generalized for vertices of higher order. However, this will not be pursued here. 



B. Trivalent vertices 



While is is not possible to decompose odd valence vertices into over and under crossings in the same manner as the 
last section, it is possible to learn about the invariants and, in particular, framing. Generically, a diffcomorphism with 
a non- vanishing volume factor will give a non-zero contribution to the variation. Take, for example, the diffcomorphism 
which rotates the trivalent vertex shown in Fig. (0). 




FIG. 4. A diffeomorphism on a trivalent vertex which corresponds to a A-move. The response of the expectation value to 
the rotation depicted by (b.) of the vertex (a.) is calculated in the text. Another embedding of the 3-valent vertex is shown in 
(c.) 



6 This may be derived by requiring that the the intersection ^ j built from ^ and , is compatible with both the 
Mandelstam identities and the Kauffman bracket skein relations (Eq. (|2^)). It might be possible to generalize this construction 
for higher valence vertices of oriented graphs. 
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One may begin with a twisted 3-vertex and rotate through an angle of 2ir. It is easy to see that the affect of this 
diffeomorphism is to cross the edges e\ and &i. Making use of the transverse 4- vertex decomposition, 




a+b 



i—\a — b\ 




(21) 



2A- 



( _ 1)( a+6-c)/2 [fl(a + 2) + b{b + 2) _ c(c + 2)] 



in which the identity of Eq. (|3J|) is used in the second line. This result is exactly what is expected from the the first 
order expansion of the A-move (Eq. (Q)! It is easy to see that the classical group recoupling determines the first 
order result. What is more, the operator is diagonal on this vertex. The eigenvalue may be then exponentiated to 
recover the whole series 



= A 



ah 



ith 



^ _ /_^\(a+&-c)/2^[a(a+2)+6(6+2)-c(c+2)]/2 



The vertex is not completely simple, though. The tangent space structure affects the result strongly. Under general 
variation the 3-vertex is 



Mi,3)+«(3,i))j>ir ( " 




9(a, c, i) 



(k(2,3) + «(3,2)) 



£(-ir 



(M 




0(6, c,i 



+ («(!, 2) + fC ( 2> l))£(-l)»'(».»>i^ 
/ — 1 V(a,b,i) 




(22) 



To reproduce the results derived above one can consider decompositions which have, for some value of u, an intersection 
between the edges e\ and e-i (and which do not move the edge 63). 

For the rotation shown in Fig. (^ (b.)) the edge remains fixed. Thus e| = 0. With the clockwise convention 
shown in the figure, this reduces to simply^] 



n in 

du— (U lu U 2n U 3 ) = - — 
du 2k 



(_i)«cK VM) + (-i)°° (6 ' c WM 



This calculation lacks the explicit projection dependence of the first calculation in that the variation is carried out 
with a diffeomorphism in the three-manifold, without reference to any preferred direction. In fact, the result of the 
diffeomorphism may be entirely different. For instance, if the first and third edges were part of a single line so that 
el = ef, as in Fig. ([|(c.)), then this diffeomorphism would only rotate the second edge. The recoupling coefficient 
would be 



2k 



k(2,1)(-1) q ^' c) &(6 + 2) 



- quite different from the previous result! This is another example of how the tangent space structure determines the 
invariant. The usual A-move is only recovered for non-collinear, essentially planar graphs with a global projection. 
(By essentially planar diagrams I mean projected graphs with intersections created in the projection labeled with 
"over" and "under" crossings.) This is expected for Temperley-Lieb recoupling theory and the Kauffman bracket as 
these are invariants for graph, knot, and link diagrams. In the next section, this becomes quite clear when single lines 
are analyzed. 



7 Here, the volume factors have been "regulated" so that each factor has the same value in the limit when the edge parameter 
for ei (and e%) vanishes as its value for any finite parameter. 
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C. Single line framing 



Framing was introduced to cure the ambiguities in the volume factor for knots. When studying invariants of spin 
nets, framing must be invariably studied on single edges. After recovering the result of Temperley-Lieb recoupling 
theory, I will examine the ambiguity of the volume term in more detail. For an embedded graph, there are several 
options, framing cycles in the graph, reaching further into tangent space, and balancing with another network as in 
the Yetter-Barret-Crane spin network invariants J33|. 

Curls may be treated as an application of the 4-vertex decomposition. In blackboard framing, since twists are 
projected as curls, the number of curls in a knot diagram is equivalent to the self-linking number. This number may 
be changed with a decomposition parameter. I take the parameter u to interpolate between curls of different chirality, 
i.e. as u flows from positive to negative, the positive curl *p changes into a negative one Jo . The intersection forms 
at u = 0. This change in the line is easily captured by the transversal 4-valent vertex decomposition. 

The variation of the parameter u results in 



-/b 




ni f i\a\^( ,^L,^ c(c + 2)\ A c 



w 



The identity of Eq. ( |45|) was used in the second line. 

There is a new feature in the calculation. To relate a curl to a uncurled line, one must shrink the loop to a point. 
One might try to perform this transformation with a second parameter much as was done for the 6-valent vertex. 
However, the volume factor vanishes. (The volume vanishes for both the first line, no intersections, and for the last 
line, a planar deformation, of Eq. (^3|).) One must merely collect the group factors. Making use of the identities of 
Eqs. pi), pi), and pi), 



(24) 



l P )-(p ) = ^(-irE(- 1 ) f |f( a ( a + 2 )- c ( c + 2 )/ 2 ) 

= (-l)*a(a + 2)j(j 

This is the expected result from Temperley-Lieb recoupling theory, ( l p ) = ( — i)°^4 a ( a + 2 ) ( j ). Again, this holds only 
for blackboard framing. 

This is a special case. One could also perform a twist on a single line. Since the boundary values of the edges must 
be left unchanged, the edge may only be rotated by even multiples of tt. Rotating one end 2ir with respect to the 
other gives 

dU ii {Ue) ^^Y I dU I dt I dse ^°(t)e 6 (i)e c ( S )<^ 

While the recoupling gives an overall factor of (— l) a; ^ra(a + 2), the volume factor is clearly ill-defined. Without 
self-intersections in the region lifted by u, there are only the one dimensional solutions s — t. Usually one regulates 
this volume factor by displacing one of the edges. By so changing one of the tangents in the volume factor of Eq. 
the delta-function may become well defined. For loops this becomes the linking number of the loop and its frame; 
to account for this framing one may associate an integer to each component of a link. In the approach taken in 
this paper, in which subgraphs with vertices are analyzed, the framing is effectively used only in a small region. As 
the self-linking of edges is meaningless, the usual approach can only be recovered if the frame is "matched" at the 
boundary. 

Alternately, one may regulate the open edge volume factor by reaching deeper into tangent space. Introducing a 
background metric and expanding one of the tangents, one finds 

(-1) — a(a + 2) / dte abc - 



k " ' 'J UUL \e{t)\ 
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which may be interpreted as an infinitesimal form of the self-linking number. 

Another point of view is suggested by the definition of a holonomy. The path ordered exponential may be given by 
the limit 



V exp 



dtd a (t) A a {a(t)) 



N 



= lim T\{l + A a da a ) 



The loop a only becomes smooth in the limit. Since it is expected that space is only effectively continuous at a 
macroscopic scale, the holonomies ought to be constructed from Planck length segments and are only defined for large 
but finite N. In this case, the effective volume factors are a set of paired terms like those given in Eq. (0). Since 
these come in pairs which only differ by sign, they cancel. In this interpretation, the difficulty of single line framing 
is an artifact of the use of a smooth manifold. It would be useful to see whether the framing difficulty remains when 
calculated on other manifolds. 

As final observation of single-line framing, it is interesting to note that the volume ambiguity may be handled in 
yet another way. Instead of regulating the delta- function and leaving the theory with an arbitrary element, one may 
instead ensure that this element cancels with another similar term. The invariant then becomes ambient isotopic 
rather than regular isotopic. To see how this may be accomplished, consider a path a and its associated frame a/. 
The framed path is related to its partner by a direction field at = a + r\9. (To capture framing information one 
is really interested in the equivalence class of smooth deformations along the loop of such direction fields. However, 
this will not affect the calculation here.) The affect is, of course, to consider two copies of the graph, one "slightly 
displaced" from the other. Typically, one frames a loop in this manner, calculates the expectation value and then 
takes the limit r\ — > 0. The result is taken to be the expectation value of the original path. Here, the variational 
method is used to see how it is possible to remove the framing dependence. 

Under a variation of a framed path both the path and its frame are parameterized by u. The variation has a total 
of four terms. 

^( U " U °t) = ^ / dt J dse abc d a ( S )d b ( S )a c (t)5 3 (a(s),a(t)) {U a (0,s)T l U^ 

dt [ dse abc d a (s)a b (s)a c f (t)S 3 (a(s), a/(t)) (U a (Q, b)T*17 (s, 1) U af (0,t)T}U af (t, 1)) 

e J a 

+ I dt [ d S e afcc d?( S )^.( S )d c (/)5 3 (^ 



a <J a 



(25) 



dt / dse o6c d?(s)d£(s)d^(i)5 3 (a^ 



Two terms have insertions in only one path, while the other two terms have insertions in both. I have omitted two 
terms like the first and last lines, which have s > t. With an eye to the limit r\ — > 0, it is reasonable to take a = d/; 
the decomposition parameter derivatives are in the same direction. This condition ensures that the second and third 
terms of Eq. ( p5| ) differ by a sign and cancel. Expanding the remaining terms, the framed path ay in terms of the 
base loop a and the frame field 9, one finds that the terms linear in 77 cancel and the last equation reduces to 

■^(U a U af ) = ^JdtJ dse abc d a ( S )d b (s)a c (t) [S 3 (a(s), a{t)) (U a (0, s)T l U a ( Sl t)T l U a {t, 1) U af ) 

+ S 3 (a f (s),a f (t))(U a U af (0,s)TjU af (s,t)TjU af (t,l))] ■ 

Despite the remaining ambiguity in the volume factor, the result of the variation can vanish if the recoupling on the 
two edges differs by a sign. This opens up a number of possibilities. One could simply insert an i in the definition of 
the holonomy. This unfortunately means that the holonomy is no longer gauge invariant in that it no longer satisfies 
the Gauss constraint of canonical quantum gravity fl4| . A far more elegant solution is to require that the variation 
of two identically labeled ST/ (2) networks differs by a sign. Since these results are first order, it suggests that one 
need only change the overall sign of the action so that e goes to e~ tS . The first order coefficient is, in the full 
series, exponentiated to the complex quantity q, so this sign change is a matter of taking the complex conjugate of the 
variable of the polynomial invariant. This is the balanced SU(2) x SU(2) invariant of Barrett-Crane p0[ and Yetter 
II- 

One might also try to build a balanced SU(2) x £7(1) network. As was noticed some time ago in the context of 
simple loops (see, for instance, J34|), the frame on a SU(2) knot can be balanced with a frame on an identical U(l) 
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knot so that the complete polynomial invariant does not depend on the frame. To see how this might arise, consider 
a SU(2) x U(l) Chern-Simons theory. The action for the U(l) part is simply 

s(a) = f d 3 xe abc a a d b a c 

in which a a is the U(l) connection. The action of the composite theory with the connection A\ = A l (T 1 )^ + aS^ 
has the same form as the SU{2) action of Eq. ([!]) p5| . Under variation the U(l) part of the theory is identical except 
for the group structure. For instance, the curl relation is 




Clearly, if the total theory is to be frame independent, the expectation values based on the two parts of the theory 
must be balanced J33|. It seems that this may be accomplished in general only if one is willing to assign irrational 
charges to the ?7(1) edges; in the general case one must require that the label on the abelian network be y / n(n + 2), 
where n is the label on the corresponding SU (2) network. 



D. Higher valence vertices 

For higher valent vertices the same recoupling formula of Eq. ( p"8| ) may be used iteratively to find the action of 
the gauge invariant operators in the spin net basis. However, there does not seem to be a canonical form of the 
intertwiners which form an eigenspace of the operators.^ Nevertheless, the recoupling procedure can be applied to 
the higher order decompositions of Eq. (fl6|). 

Despite this, there is a class of graphs which have a particularly simple decomposition. These graphs have edges 
in the fundamental representation (labeled by 1) and even-valence vertices which are "consistently oriented." That 
is, all the vertices have incident edges which are oriented in alternating directions as projected along the vector field 
associated which the decomposition (as in Fig @). For these graphs, the evaluation of the operators, making use of 
the identity of Eq. (^), yields simple connection diagrams; the variation relates the decomposition of the vertex into 
over and under crossings to a sum of links with non-intersecting components. Assuming it is possible to "smooth" 
the edges to remove the kinks at the vertices, the affect of the operators then simply produce a sun of regular link 
invariants. In these cases, the Vassiliev invariant is given by the chromatic evaluation of the loops. 



V. DISCUSSION 



In this paper new invariants for embedded graphs in Chern-Simons theory. The key difference from earlier work 
is that the new invariants depend on the tangent space structure at vertices. Using the variational technique, this 
analysis of the definition of the vacuum expectation value of embedded graphs or, equivalently, of the spin net 
representation of the Kodama state, suggests that one may sensibly define invariants of graphs. For closed graphs, 
these operators are simply a sum of signed terms consisting of Wilson graphs with generators inserted at the vertices. 
The result is a Vassiliev invariant of order (n — 1) for a 2n- valence vertex and so provides an example of the theorem of 
Birman and Lin on the expansion of a polynomial invariant. The beauty of the variational method of graph invariants 
lies in that it creates relations between different non-intersecting invariants of the 3-manifolds without resorting to 
a fixed background structure. The result is an invariant which depends not only on the 3-manifold but also the 
tangent space at the vertices. Graph invariants then capture information of the three manifold as well as the tangent 
space at the vertices. In the spin net basis it is possible to evaluate the action of these operators. Some of these 
variational operators may be formally exponentiated to give the result to all orders. In this manner, the Temperley- 
Lieb recoupling theory of Kauffman and Lins is recovered for essentially planar diagrams with blackboard framing. 
The calculation reveals that it is not necessary to separately frame vertices as the variation is well-defined without 



If the recoupling identity used in calculating the area operator eigenvalues in the diagrammatic approach could be generalized 
for the coefficients a c (a,b), it would easy to find a suitable intertwiner tree. Alternately, there may be a relation among the 
spin operators at the vertex which suggests an intertwiner basis. (See Refs. 
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further regulation. By examining from this perspective the fra ming of a single edge of a spin network, the balanced 
spin network invariants of Barrett-Crane- Yetter were recovered f20f , [j33f . 



Motivated by considerations arising from canonical quantum gravity, this study focused on the role of framing at the 
vertices. Though spin net geometry does not require framing to be rigorously well defined, there are three immediate 
reasons why one might expect that framing is a property required by the full theory. First, the cosmological constant 
appears in the invariant. Since the cosmological constant appears only in the Hamiltonian constraint, framing is an 
issue of dynamics. While it may be that framing is only required for the "Kodama phase" of the theory, the state 
suggests that the complete theory (taking seriously the suggestion that the Kodama state is a well-defined solution 
of the full theory) will need to account for framing in at least one sector. Second, in the loop representation of the 
linear theory |36| as well as Maxwell theory [p7| , framing plays a key role. Third, since framed links are sufficient to 
construct all compact, oriented 3-manifolds |38[ , a theory of the dynamics of such manifolds ought to have a framed 
loop representation. 

Even so, the Kodama state provides an enigma for canonical quantum gravity. While it has all the expected 
characteristics and is a formal solution to the constraints, it cannot be treated with the theory of measures which 
has proven so fruitful in spin net geometry . A key conceptual confusion has been the lack of understanding of 
the requirement of diffeomorphism invariance. Is the classical three-manifold diffeomorphism invariance broken by 
dynamics? In the language of knot theory, are kinematic states regular or ambient isotopy invariants? While it seems 
obvious that when dealing with such a canonically diffeomorphism invariant theory as gravity that we must consider 
invariants of ambient isotopy, such as the Jones Polynomial, it is only clear that in the classical limit gravity possesses 
the full diffeomorphism invariance. Indeed, given the universal nature of Chern-Simons action, it even seems possible 
that this state may be the scaling limit of an underlying theory of quantum gravity. It perhaps indicates that a more 
sensitive invariant may be required to describe the full theory. One can only hope that it is pointing to a feature of 
the microstructure of spacetime. 

On a more prosaic level, the representation of framed spin networks created in part to describe states in the 
Kodama phase pH , offers a number of simple mathematical challenges: What is the product of two framed loops 
which intersect? Is such a product unique? What is the appropriate algebra for the basic operators of the theory? 
The variational and recoupling techniques offer a method for resolving these issues. For these reasons, this study 
concentrated on the issues of frame and decomposition of embedded graphs. 

The results suggest that is is possible to further and consistently define the vacuum expectation value of a graph. 
Through a delicate interplay of group and tangent space structure, the invariant seems to be defined. It is reasonable 
to expect that framed graphs include tangent space information in addition to that of the three manifold. As is 
suggested by geometric operators in spin network kinematics, the tangent space plays a large role. The invariants of 
framed graph depend critically on the tangents of the incident edges. This work offers a hint of how tangent space 
and manifold information might be folded together. One potentially interesting direction to explore is the evaluation 
for flat connection. Since the holonomies reduce to identity, the result would be a numerical knot invariant p2[ . 
Composed of group, manifold, and tangent space structure, these invariants would give some intuition of spin net 
expectation values as well as the state space of spin net geometry. 

I would like to close with two additional comments. The first is a bit technical. As was noted in the Section III, the 
form of the vertex decomposition is given by the choice of the decomposition vectors at the vertices e a \ u=Ua . These 
vectors determine both the signs arising from the volume factor and which non-intersecting links appear. The freedom 
in the choice of these vectors is the freedom in the definition of the vertex |3^] . One way to identify this freedom is by 
identifying the "tangent space volume" in which the decomposition vector lies |52j . This is based on the observation 
that a diffeomorphism acts as an invertible linear transformation on the incident edge tangents. A volume, identified 
by a non-co-planar triple, cannot be made to vanish by a diffeomorphism. This suggests that the amount of freedom 
in defining a vertex is contained in how the graph is embedded in the manifold. 

The second comment is more speculative. Since there is no natural relation between the decomposition vectors at 
distinct vertices, the graph invariant seems to be ambiguously defined without a global projection. If the expectation 
values of spin nets are invariants of the 3-manifold rather than invariants of planar diagrams, then there must be a 
consistent definition. As the expectation value may differ by framing and as the framing dependence collects into an 
overall factor, I would like to suggest a gauge principle for the framing of embedded graph invariants;^ it ought to 
be possible to choose any decomposition of each distinct vertex in a graph and still find the same invariant - up to 
"gauge." There is a hint of this in the planar representation of framed graphs. One simple example is of two subgraphs 



9 This idea was first mentioned in a discussion with Lee Smolin, Roumen Borissov and myself at the Center for Geometry and 
Gravitational Physics. 
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connected by a single edge. In a global projection, distinct decompositions of the graph ought to be identical up to 
gauge. For instance if one subgraph is rotated 2n with respect to the other, the result is, of course, a curl. This 
may be the phase of an abelian gauge theory. Another simple example is given by a pair of trivalent vertices joined 
by a pair of edges. In a global projection, the decomposition of the graph ought to be identical, up to gauge, if one 
vertex is rotated 2ir about an axis joining the two vertices. A simple application of recoupling theory shows that 
the affect of the rotation is simply equal to a curl on one of the "external lines." This again suggests a curl as a 
"phase" due to a gauge rotation. While it is not clear from these examples whether such a gauge principle would 
hold on arbitrary vertices and arbitrary decomposition vectors, it nonetheless suggests that the framing dependence 
of vacuum expectation values of embedded graphs in Chern-Simons theory is gauge. 
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VII. APPENDIX 



This appendix contains the basic definitions and formula of recoupling theory. While this work uses for the most part 
the binor conventions for SU(2) recoupling theory, the formula here are for more general A (except where explicitly 
stated). For more than a brief review see Ref. JlSj] . The complex phase A is given by 

A = e l7I,2k 

for integer k. A is found in the fundamental skein relation for the Kauffman bracket 

X =^X ■ ( 2? ) 

and is related to the usual parameter q via 

q = A 2 = e™' k . 

In, g-quantum gravity, this parameter is given by 




so that k = 6n/Alp (an integer!). F] The "classical" or "binor" limit occurs when q — 1 and A = — 1 (r — > oo) so that 
h and/or A vanish. In this limit, the relation of Eq. (^) (which is also the relation among 2x2 matrices known as 
the Mandelstam identities) take the simple diagrammatic form of 

X + )( + X =°- ( 28 ) 

Recoupling theory begins with the basic irreducible representation, diagrammatically a single line or "strand." 
Closing this line (or taking its trace) gives the loop value d = — A 2 — A~ 2 . The classical value (when A = ±1) of d is 
—2. Higher representations may be built from the basic line using the Wenzel- Jones projector defined by 

A =tXE(^ 3 ) H ^ (29) 

| 1 ; '<tGS„ 



10 If one includes CP-breaking term in the action, f F A F, and the non-perturbative renormalization parameter jj, then one 
has k — > k + 2 with 

, 6-7T 



The parameter a is the phase coming from the CP-breaking term. 
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in which the sum is over elements of the symmetric group, er; \a\ is the sign of permutation; the expansion c-n is 
given in terms of the positive braid (the strands are only over crossed )\ ); and the asymmetric quantum number 
{n} is defined by 



{n}:= 



1-A- 
1-A 

This quantum integer is simply an integer in the classical limit. 
The evaluation of a single un-knotted n loop is 



-4 ' 



(-!)>+!] 



(30) 



(31) 



where [n + 1] is the dimension of the representation and the brackets identify the symmetric quantum integer defined 
by 



A 2n _ A 



-2n 



11 A 2 -A- 2 ' 

For classical spin networks A = — 1, the next two identities are useful in the calculation of the curl 

2a 

£ (-l)«A c = A a 



(32) 



and 



til! I <_ 1)S £<£±^ = a(a + 2 , 

c— 0;cvcn 



The function 9(a,b,n) is given by 

9{m, n, I) 



(-1) 



(a+b+c) 



[a + 6 + c + l]![a]![b]![c]! 
[a + b]\[b + c]\[a + c}\ 



(33) 



(34) 



where a = (I + m — n)/2, b = (m + n — l)/2, and c = (n + I — m)/2. A "bubble" diagram is proportional to a single 
edge. 



(-i)"fl(a,M) 
[n + l] 



(35) 



The basic relation relates the different ways in which three angular momenta, say a, 6, and c, can couple to form a 
fourth one, d. The two possible recouplings are related by the formula: 



c 



E 



\a-b\<i<(a+b) 



a o i \ ! s r\ 
c d i' | \ X. , ! 



(36) 



where on the right hand side is the (76j-symbol is defined below. It is closely related to the Tet symbol. Variously, 
drawn as 



it is defined by ml 



bee 



b c e [,' 

or /-> ? 

d 



- f J} I = Tet 

d 



a b e 
c d f 



Tet 



a b e 
c d f 



n y (-iY iiril] 



(37) 



N 



m<B<S 

[a]\[b]\[c]\[d]\[e}\[f]\ 
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in which 



ai = ^(a + d + e) 
ai = \{b + c + e) 
a 3 = i(a + 6 + /) 
a 4 = §(c + d + /) 

The <?6j-symbol is then defined as 



b 1 = \{b + d + e + /) 

62 = i(a + c+e + /) 

63 = i(a + 6 + c + «i) 
m = max {a^} M — min {bj} 



a b i 
c d j 



These satisfy a number of properties including the orthogonal identity 



Tet 


a b i 
c d j 


A, 


9(a, c 


l,i) 9(b, 


c,i) 



E 



a b I 
c d j 



d a i 
b c I 



and the Biedenharn-Elliot or Pentagon identity 

E 



d i I 

e m c 



a b I 

e I i 



a f k 
d d I 



a b k 
c d i 



k b f 

e m c 



The "A-move" 



A. 



A" 



where A" fc is 



^ab _ (_ 1 )(a+&-c)/2 j 4[a(a+2)+b(6+2)-c(c+2)]/2 

An over-crossing may be related to a recoupling via 



a+b 



= E ^ 



6*(a, 6, i) 



i— |a — 6 

A similar identity holds for undcr-crossings. In the binor limit the two recouplings coincide. 



(38) 
(39) 
(40) 



(41) 



(42) 



(43) 



(44) 



(45) 
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